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Abstract 



The symmetry-induced magnetic anisotropy due to monoatomic steps at 
strained Ni films is determined using results of first - principles relativistic 
full-potential linearized augmented plane wave (FLAPW) calculations and 
an analogy with the Neel model. We show that there is a magnetoelastic 
anisotropy contribution to the uniaxial magnetic anisotropy energy in the 
vicinal plane of a stepped surface. In addition to the known spin-direction re- 
orientation transition at a flat Ni/Cu(001) surface, we propose a spin-direction 
reorientation transition in the vicinal plane for a stepped Ni/Cu surface due 
to the magnetoelastic anisotropy. We show that with an increase of Ni film 
thickness, the magnetization in the vicinal plane turns perpendicular to the 
step edge at a critical thickness calculated to be in the range of 16-24 Ni layers 
for the Ni/Cu(l,l,13) stepped surface. 
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Ni/Cu(001) is known as a unique system to show an in-plane to out-of-plane spin- 
reorientation transition with an increase of Ni-overlayer thickness |jl|J^. This is hkely due to 
the positive magnetoelastic energy contribution to the magnetic anisotropy energy (MAE), 
caused by tetragonal strain in Ni- films due to the film/substrate lattice mismatch 

We have investigated the uniaxial MAE for flat and stepped Ni films, and use the results 
of first-principles MAE calculations for the strained Ni-films to estimate the in-plane to 
out-of-plane spin reorientation transition thickness. Further, we also use an analogy with 
the Neel model with parameters chosen from our first-principles calculations to estimate the 
step-induced contribution to the MAE and to discuss the possibility of a spin-reorientation 
transition at the Ni/Cu(001) monoatomic step. 

The relativistic self-consistent version of the full-potential linearized augmented plane- 
wave (FLAPW) method is employed to perform self-consistent calculations for strained 
Ni films and bulk Ni with an (001) magnetization direction. The local force theorem for 
magnetization rotation is then used for the MAE calculations [0]. 

For strained bulk Ni, we choose the bet (001) plane lattice constant as 4.781 a.u. (equal 
to the Cu substrate) and the (001) interlayer distance as 3.222 a.u. in accordance with the 
experimental strained structure of a Ni magnetic overlayer on Cu(OOl) ||^. Here, 343 k-points 
in l/8th of the tetragonal 3D BZ were used for the self-consistent calculations. Then 1728 
k-points in l/8th of the 3D BZ (13824 in the full BZ) are used for the MAE calculations. The 
MAE is shown in Table |. There is a very good agreement between our results and results of 
previous perturbative state-tracking FLAPW calculations and relativistic full-potential 
linear muffin-tin orbital (FLMTO) calculations [Q. 

For Ni-films with a thickness (d) of 3,5,7,9,11 layers we choose a slab geometry with 
assumed experimental strained structure of Ni magnetic overlayer on Cu(OOl) which is the 
same as in the strained bulk Ni calculations. Here, 100 k-points in l/4th of the 2D BZ (400 
k-points in the full 2D BZ) were used for the self-consistent calculations and 1600 k-points 
(6400 k-points in the FBZ) in 1/4 of the 2D BZ were used for the MAE calculations. 

The calculated MAE as a function of Ni-film thickness is shown in Fig.l. There is 
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a very good linear dependence of the MAE on film thickness. The uniaxial MAE of d 
magnetic bct(OOl) layers in the [(100), (010), (001)] coordinate frame (denoted by V) is then 
approximated by 

E =-d*K^*V^ -2*Ks*V^ (1) 

where, is a "volume" MAE per atom. Kg is a "surface" MAE per atom, and Vz is 
the magnetization direction cosine with respect to the z(V)-axis. By fitting the calculated 
MAE (cf.. Fig. 1) by Eq.(|l]), we obtain the "volume" = 83.5 fieV and "surface" 
Kg = —446.5 fieV contributions. The calculated "volume" contribution is in very good 
quantitative agreement with the calculated uniaxial MAE for strained bulk-Ni (cf.. Table |). 
This is clear numerical evidence that the "volume" contribution to the MAE in strained Ni 
films is the same as the MAE of strained bulk Ni, which is qualitatively consistent with Neel 
model arguments [Il5| , [l6| and supports quantitatively an analysis performed in Refs. |]l],§]. 



Hjortstram et al. 0] have calculated the uniaxial MAE and magnetostriction in strained 
bulk Ni and found good agreement between the linear magnetoelastic theory and the results 
of first-principles calculations for strain-induced MAE in the limit of small strains. They 
also found that for assumed experimental strained structure of Ni films the MAE value for 
strained bulk Ni agrees well with measured "volume" contribution to the MAE for Ni films 
on Cu(OOl) ||l|]. The present calculations for both strained bulk Ni and Ni films (cf.. Table 
I) agree well with the results of The fact, that "volume" MAE (85 ^eV) is somewhat 
bigger than the MAE for strained bulk Ni (57 ^eV) can be due to computational differences 
and neglect of the higher order terms in MAE expression Eq.(l) (which are small but not 
necessarily negligible). Taking into account relatively small calculated MAE values we argue 
that both the "volume" MAE for Ni films and MAE for strained Ni bulk agree surprisingly 
well. This good numerical agreement makes it possible to conclude that "volume" MAE 
in Ni films is due to the strain induced by film/substrate lattice mismatch. This confirms 
quantitatively conclusions of Hjortstram et al. [Q based on comparison between experimental 
data for fiat Ni films and the results of strained bulk Ni calculations. 



Finally, using Eq.([l|) and the calculated MAE data (cf., Table HI), we calculate the in- 
plane to out-of-plane spin-reorientation transition thickness as dc = 11 ML, which agrees 
qualitatively with the experimentally observed spin-reorientation transition thickness of 7 
ML for Ni/Cu(001). (The demagnetization energy is omitted here because it is very small 
for Ni 0.) It is necessary to mention that the use of free-standing Ni- films to analyze the 
MAE for Ni/Cu(001) is rather qualitative since we do not take into account the interface 
contribution to the MAE and the possible oscillations of the MAE with Ni-film thickness at 
the Ni/Cu(001) interface Moreover, the chosen geometry is an "average" experimental 
geometry and deviations due to relaxation are likely to occur in the real interface. 

To proceed further, let us now assume the Neel pair - interaction form of the mag- 
netic anisotropy energy {Eab) for the pair of atoms A and B connected by Hab = 
I'^abUPx, Py, Pz), with magnetic moment direction m = (yx,Vy,Vz), where P,V are the 
direction cosines; Lab is a Neel coefficient (the isotropic term is omitted). Hence, 

Eab = - Lab{ E AV^.)' (2) 

i=x,y,z 

The "volume" and "surface" MAE contributions in Eq.(|]) are then expressed as: 

K = +L{Q(3l - 3) (3) 

Ks = -L{3P! - 1) (4) 
where, following Victora et al. we consider only nearest-neighbor interactions and do 



not take into account the change of the Neel parameter due to the strain (which seems to 
be reasonable assumption in the limit of small strains, considered here). 

The value of the Neel parameter, L = —396 fieV, is then calculated using Eq.(|), the 
calculated value of the "volume" anisotropy for strained bulk-Ni (cf.. Table D and = 0.476. 
The value of the "surface" anisotropy is then calculated using Eq.(^) and is 169.5 fieV. A 
comparison with the result of first-principle calculations (cf.. Table D clearly shows that this 
"isotropic" variant of the Neel model does not even reproduce the sign calculated from the 
first-principles "surface" MAE. 



It is known that there is a significant change in the electronic structure at the Ni- 
surface compared to the bulk. This change is especially pronounced for the "surface" layer 
and leads to significant modifications of the electronic density of states near the Fermi level 
for the spin-minority band. This should lead to a significant difference in "in-plane" and 
"inter-plane" Neel parameters for the surface layer. Therefore, it is quite natural to introduce 
an additional in-plane Neel parameter M to describe the pair-interaction, Eq.(0), within the 
surface layer. This gives the following parametrization for the "surface" anisotropy: 

= -L{3(3l - 1) + (L-M) (5) 

Then, we approximate the MAE of a strained Ni-film by choosing the Neel parameter 
L to fit the "volume" contribution to the MAE, Eq.(^, and then determine the in-plane 
parameter M to fit the "surface" contribution, Eq.(|^). Using the calculated values of the 
"volume" and "surface" anisotropics for the strained Ni-film (cf. Table ||), we obtain the 
following values for the Neel parameters, L = —580 fieV and M = 115 fieV. 

As shown in Refs. [0,0], the anisotropy due to the step along (110) on the fct(OOl) 



surface can be written in terms of the contributions from "step - corner" and "step - edge" 
atoms (we use the bet-system and a fct (110) step is a bet (100) step): 

Est = -KseV^ - KlV^ - KlV.V. (6) 

where, Kse is a "step-edge" atom's MAE, and K]^ and K^^ are a "step-corner" atom's MAE. 

Within the assumption of only nearest-neighbor interactions one can express i^se? 
and K^^ in terms of the Neel parameters L and M: 

Kl = 0.51(3^', - 1) 
Kl = -LM. 

Kse = -0.5M (7) 

The values of Kse,Kl^ and K]^ are then calculated for the Ni- step using Eq.(^ and are 
shown in Table |I| together with values of the "volume" and "surface" anisotropics. 



The symmetry induced uniaxial one-step MAE (per one atom in the xy-plane) is then 
nven 
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E = -D * dK^V^ -2* D * KsV^ - KseV^ (8) 

where D is the number of one-step atoms (along x). We consider only the "surface" step 
and neglect the "interface" step. As follows from Eq.(^, the direction of the magnetization 
in the xy-plane is determined by the sign of Kse- In the case of a Ni stepped surface, 
the magnetization is directed along the step-edge (perpendicular to the x-axis) due to the 
negative sign of Kse (cf. Table This is in qualitative agreement with experimental data 
for thin Ni-films on stepped Cu(OOl). 

Next, we transform Eq.(|^) to a vicinal coordinate frame (U) using the relation between 
the old (V) and new (U) magnetization direction cosines in terms of the vicinal angle a. 
Keeping terms up to a^, and with 

K = (1 - 0.5a2)f/^ + aU, (9) 
v; = (1 - Q.W)U, - aU^ 

we come to the following expression for the MAE in the xy-vicinal plane 

= {-D*dK^a^ -2D*Ksa^ -Kse{l-a^) (10) 



Using the relation between a and D (a ~ lop/{aipD) = x/D) and neglecting 1/D^-order 
terms we transform Eq. ([lO|) into: 



E'iy = i-dx^D K,-2xyD Ks-Kse (11) 



-x/DKl) Ul 



We next apply Eq. ([TT|), with the parameters determined within the framework of the 
modified Neel model (cf.. Table |I|), to calculate the MAE at the (1,1,13) stepped surface. 



Using the strain parameter x = 0.674 and step lenght D = Q, we obtain the following 
contributions to the MAE: 

(1) The contribution due to the "volume" magnetoelastic anisotropy is —6.322 x d fieV; 
it turns the magnetization perpendicular to the step and increases with an increase of the 
film thickness; 

(2) The contribution due to the "surface" MAE is 68 fieV; it turns the magnetization 
parallel to the step and does not depend on the film thickness; 

(3) The contribution due to the "step-edge" MAE is 57.5 fieV; it turns magnetization 
parallel to the step and does not depend on the film thickness; 

(4) The contribution due to "step-corner" MAE is +22 fieV, and it turns the magneti- 
zation parallel to the step and does not depend on the film thickness; 

Finally, for thicknesses below 24 ML, the contributions (2), (3) and (4) will keep the 
magnetization parallel to the step in the vicinal (xy)-plane and the critical thickness for the 
magnetization to turn perpendicular to the step edge is then estimated to be 24 ML. 

In conclusion, Eq. (|ll|) shows that there is a "volume" uniaxial magnetoelastic anisotropy 
contribution to the symmetry-induced magnetic anisotropy in the vicinal plane of the 
stepped surface, and it is caused by tetragonal strain in the magnetic film due to 
film/substrate lattice mismatch. For monoatomic steps on Ni-films, this magnetoelastic 
contribution leads to an additional spin-reorientation transition in which the magnetization 
turns from parallel to perpendicular to the step alignment with increase of the Ni-film thick- 
ness. Using the results of first-principles calculations for the strained Ni film and a Neel 
model-type analysis, we estimate the critical Ni-film thickness for this transition to be 24 
ML for the (1,1,13) stepped surface. One has to take into account the difference between 
the surface anisotropy in Ni/Cu(001) films and in the calculated free-standing Ni-films: if 
we use the experimentally derived (at zero temperature) "surface" MAE ^ (i.e. -250 fJ^eV), 
we expect a critical thickness of 16 ML for this transition. This result shows the uncertainty 
and range of values possible at this level of theory. 
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TABLES 

TABLE L Uniaxial MAE of tetragonal strained bulk Ni (K^), "volume" (K^) and "surface' 
{Ks) contributions to the MAE for strained Ni-films. All data are given in ^eV. 
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TABLE IL Values of K^, Kg, Kse,Kl^ and K"^^ (in ^eV) for the Ni-step. 







Kse 




sc 


+83.5 


-446.5 


-57.5 


-124 


195 



11 



FIGURES 

Fig. 1 The MAE (in n eV) for the strained Ni-films as a function of the film thickness 
d (in ML). 
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